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Abstract 

We study tachyonic preheating associated with the spontaneous breaking of B — L, 
the difference of baryon and lepton number. Reheating occurs through the decays 
of heavy Majorana neutrinos which are produced during preheating and in decays of 
the Higgs particles of B—L breaking. Baryogenesis is an interplay of nonthermal and 
thermal leptogenesis, accompanied by thermally produced gravitino dark matter. 
The proposed mechanism simultaneously explains the generation of matter and dark 
matter, thereby relating the absolute neutrino mass scale to the gravitino mass. 



1 Introduction 

Thermal leptogenesis [T] explains the observed matter-antimatter asymmetry of the uni- 
verse in terms of neutrino masses that are consistent with neutrino oscillation experiments. 
In its simplest version, the primordial lepton asymmetry is generated by the CP violating 
interactions of N\, the lightest of the heavy Majorana neutrinos, the seesaw partners of the 
ordinary neutrinos. Typical values for the N\ mass and for the leptogenesis temperature 



In supersymmetric theories such high reheating temperatures cause the 'gravitino 
problem' for heavy unstable gravitinos [3|T}5] . However, if the gravitino is the light- 
est superparticle (LSP), this problem can become a virtue since for typical leptogenesis 
temperatures and for superparticle masses of the electroweak scale, thermal production of 
gravitinos can explain the observed amount of dark matter [6] . The dominant contribution 
from QCD processes is given by 



where Tr is the reheating temperature while tuq and m- g are gravitino and gluino masses, 
respectively; the coefficient C = 0.27 to leading order in the gauge coupling 

[38] Hence, 

for Tr ~ Tl and superparticle masses characteristic of gravity or gaugino mediation one 
indeed obtains the observed amount of dark matter, Odm/j 2 — 0.11 [10J. 

Thermal leptogenesis is consistent with a reheating temperature Tr much larger than 
Tl, which is not the case for thermally produced gravitino dark matter. Consistency of 
these two mechanisms for the origin of matter and dark matter then requires an explana- 
tion of why Tl and Tr have the same order of magnitude. 

Successful leptogenesis, independent of initial conditions, favours the mass window 
10~ 3 eV < rrii < 0.1 eV [H] for the light neutrino masses . For an effective neutrino mass 
rh\ = (m) D mD)ii/ Mi ~ 0.01 eV, where mp is the Dirac neutrino mass, the decay width 
of the heavy Majorana neutrino Ni becomes 



where we have used t> EW = 174 GeV for the vacuum expectation value of electroweak 
symmetry breaking. If the hot phase of the early universe is initiated by the decays of 

1 Note that C has an 0(1) uncertainty due to unknown higher order contributions and nonperturbative 
effects [7] • Resummation of thermal masses increases C by about a factor of two [5] . 



are M 1 ~ T L ~ 10 10 GeV [2]. 
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particles with decay width T, the reheating temperature is given by 




where g* ~ 200 is the effective number of relativistic degrees of freedom and Mp = 
1.22 x 10 19 GeV is the Planck mass. If reheating occurs through the decay of Ni neutrinos, 
one obtains 

T R ~ 0.2 • ^T° Ni M P ~ 10 10 GeV , (4) 

which is indeed the temperature needed for thermal leptogenesis and gravitino dark mat- 
ter. 

This observation raises the question whether in the thermal bath produced by heavy 
Majorana neutrino decays both, the matter-antimatter asymmetry and gravitino dark 
matter, can be simultaneously produced. In the following we shall present an example 
which demonstrates that this is indeed possible. 



2 Flavour model and tachyonic preheating 

In supersymmetric extensions of the Standard Model the superpotential for matter fields 
including right-handed neutrinos reads 

W M = /^KUO^ + hffilOjHi + h'^n'jll,, + Kn^S , (5) 

where we have used 577(5) notation for ordinary quarks and leptons, i — 1, 2, 3 counts the 
generations and n c denote the charge conjugate of right-handed neutrinos; H u , H d and S 
are the electroweak and B — L symmetry breaking fields, respectively, 

(H u ) = v u , (H d ) = v d , (S) = v B _ L . (6) 

The heavy Majorana neutrinos are given by Ni = rij + n\. 

For the Yukawa couplings we shall use a representative model [12] with a Froggatt- 
Nielsen U(l) flavour symmetry, which is consistent with the observed hierarchical quark 
and lepton masses, thermal leptogenesis [12] as well as constraints from flavour changing 
processes [13]. Their order of magnitude is given by 

hij oc , (7) 

where the chiral U(l) charges are listed in Table dj the parameter rj ~ 1/V300 is deter- 
mined by the hierarchy of quark and lepton masses. 
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Table 1: Chiral 17(1) charges. 10 = (g,u c ,e c ), 5* = (d°,Z). 



Using the seesaw mass relation and the estimate m u = v /m 3 m 2 ~ 0.01 eV, successful 
thermal leptogenesis determines the mass of the lightest heavy Majorana neutrino, M% = 
10 10 GeV [12] . From 

Mi = klVB-L ^ rfvB-L (8) 

one then obtains vb-l — 3 x 10 12 GeV. The masses of the other heavy neutrinos coincide 
with the B — L breaking scale, 

M 2 ~ M 3 ~ v B ^ L ~ 3 x 10 12 GeV . (9) 

For simplicity, we shall only consider the Higgs boson of B — L symmetry breaking 
in the following 2 ]. The Higgs potential relates the S mass to the energy density of the 
unbroken phase, 

m 2 s = \ V 2 B _ L , Po = ^Xv B _ L . (10) 

We assume the Higgs coupling A to be 0(1), such that the Higgs mass lies in the range 

2M 2 , 3 >m s > 2Mi . (11) 

The U(l) flavour symmetry also determines the decay rates of the heavy Majorana 
neutrinos Ni, N 2 ^ and of the Higgs boson S, 

r° Nl ^ |U , (12) 
,2 



^3^^2,3, (13) 



S7T 

= [1 - (2M 1 /m 5 ) 2 ] 3/2 ~ ^-m s . ( I 1 ) 

167T l07T 

Note that these estimates for the decay widths have a considerable uncertainty since the 
Yukawa couplings, which enter quadratically, are only known up to factors 0(1) [12] . 



The supersymmetric Higgs mechanism for B — L breaking will be discussed in a forthcoming publi- 
cation rm. 



For instance, since M\ and T° Ni involve different Yukawa couplings, the ratio T° N JM\ 
can easily vary by two orders of magnitude. Since this ratio will be important in the 
following phenomenological analysis, we shall parametrize this uncertainty explicitly by 
the effective neutrino mass fh\ introduced in Eq. (j2J). 

A crucial ingredient of leptogenesis are the CP asymmetries in heavy Majorana neu- 
trinos decays [T5fT6] . 

e< = — „ \ , Vim \(h^h v )U F (rfp) , (15) 

from which we obtain 

d ~ 0.1 if ~ 10~ 6 , e 2>3 ~ 0.1 rj 2 ~ 3 x 10~ 4 . (16) 

These asymmetries determine the lepton asymmetry which is generated by decays and 
inverse decays of the heavy neutrinos. Note that for Mi ~ 10 10 GeV, the estimate 
ei ~ 10~ 6 corresponds to the maximal possible CP asymmetry [17J. For smaller heavy 
neutrino masses the asymmetry is reduced by the mass ratio, e% ~ 10 -6 Mi/10 10 GeV. 

The vacuum energy of a phase with unbroken B — L symmetry can drive a period of 
inflation, which then rapidly ends due to the spinodal growth of long wave-length modes, 
a process known as 'tachyonic preheating' [18J. The true vacuum is reached at a time tpn 
after inflation |19j . 

( s1s )L tm = v l-L, * PH ^ln(^). (17) 

During this phase transition the energy of the false vacuum is converted mostly into a 
nonrelativistic gas of S bosons, with an admixture of heavy neutrinos^ whose contribution 
to the energy density is determined by their coupling to the Higgs field [T9] , 

r Ni = ^ ~ 1.5 x 10- 3 g N Xf {a h 0.8) ; (18) 
Po 

here g N = 2 and / (a, 7) = a/« 2 + 7 2 — 7 with aj = hf/VX. For the neutrinos A^,3 one 
obtains 

r N2>3 ~ 1 x 10- 3 (h%) 2 ~ 1 x 10" 3 . (19) 

The fraction = 0{r] i ) and it is therefore negligible. 

3 For simplicity, we neglect the contribution of scalar neutrinos ; also in thermal leptogenesis their 
effect on the final baryon asymmetry is known to be small [20 . 
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3 Reheating, baryon asymmetry and gravitino dark 
matter from heavy neutrino decays 



We are now ready to study the time evolution of the initial state produced by tachyonic 
preheating. First, the dominant contribution to the energy density is the nonrelativistic 
gas of S bosons (cf. Eq. (Fl9l)). 



Ps{t) 



a 3 (t) 



2rjv 2 )po 



(20) 



where a(t) is the scale factor of the Friedman universe which is normalized by setting 
&(£ph) — 1- It is convenient to introduce the comoving number densities for all particle 
species X, 



N x (t) = a 3 (t)n x (t) . 



(21) 



At t ~ t 2 = 1/T% 2 the heavy neutrinos N 2 and N 3 decay, producing the initial 
radiation densit\n 



2r N2 p 
a 3 (h) 



(22) 



with temperature 



T(t 2 ) 



30 



1/4 



7r 2 g* a 3 (t 2 ) J 



(23) 



The out-of-equilibrium decay of N 2 and A^3 also produces a B — L asymmetryo The 
corresponding comoving number density is given by 



N. 



B-L 



(t 2 ) = e 2 N N2 (t 2 ) + e 3 N N3 (t 2 ) ~ 0.2 r) 2 N N2 {t 



PH; 



(24) 



The scale factor a(t) has to be determined by solving the Friedman equation. For a 
flat universe and constant equation of state uj = p/p between some initial time t and 
time t, one has 



a(t) = a(to) 



l + -(l + o;) 



S7T 



3M|" 



1/2 

Ptot (t ) ) (t - t ) 



3(1+^) 



(25) 



Since ti is very small compared to the time scales relevant for leptogenesis, we approximate the onset 
of radiation by a step function. 

5 In [19] only the initial B — L asymmetry from tachyonic preheating is taken into account. 
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One expects that until the decay of the S bosons at t ~ ts = l/^s the s yst em is mostly 
nonrelativistic and that it becomes relativistic at later times. We have checked numerically 
that this is indeed the case with effective equation-of-state parameters 3(1 + wph) — 3.3 
for tpH < t < ts and 3(1 + us) — 4.0 for t > ts, respectively. 

The thermal part of the plasma due to N 2 and N 3 decays produces also Aq neutrinos. 
Their comoving number density N~[ satisfies the familiar Boltzmann equation 



(1 rp 

aH Ta N «> 



-r« (iVj, - Aft) ; 
here H = a/a is the Hubble parameter, and the thermal width 



(26) 



(27) 



K 2 {z) ' T ' 

where K\ and K 2 are modified Bessel functions. Note that in deriving Eq. ( 126]) one 
assumes kinetic equilibrium for the thermally produced N\ neutrinos. 

Ni neutrinos are also produced in S decays. These neutrinos are relativistic, but not 
in kinetic equilibrium. They are produced with energy ms/2 which is then redshifted with 
increasing scale factor. We take this into account by solving the Boltzmann equations for 
the S and Aq distribution functions (.EWi = \lv 2 + 



--Hp- 

dt dp 

--Hp- 

dt dp 



fs(t,p) 



E s 
Mi 



E 



Ni 



nfs(t, P ) , 



(2f 



+ 



2Tt 2 n s T\ 
E 2 



[1 - (2M 1 /m s ) 2 ] 1/2 5 (E Nl - ms/2) 



A straightforward calculation yields for the distribution functions 



fsfap) 



k 2 

i 47r 2 r° 

a 3 ins 



a{t) 

t 

(l-(2M 1 /m s ) 2 ) 4 / dt' 



(29) 



(30) 



t2 



-5{E Nl -£ Nl {t\t)) 



x e-\ (t\ t) N s (f) exp f -M^ J dt"S N l (f , t") 



(31) 



where S^it 1 ,t) is the redshifted energy at time t of an Aq neutrino produced in S decay 
at time t' with energy ms/2, 

1 



£ Nl (t,t)-—- 



-<<r 



m s 



(32) 
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The corresponding comoving number densities are easily obtained from 

N x (t) = a 3 g x J -0-f x (t,p) , (33) 

where X = S,N l . 

The number densities of thermally and nonthermally produced Ni neutrinos enter as 
source terms in the Boltzmann equation for the B — L asymmetry, which reads 

7 TV 6 " 1 ~ 

*H-N B _ L = ei T Nl {Nl -N e £)- -^T Nl N B „ L + e 1 T%N s Ni , (34) 
with the nonthermal contribution 

Nm-)dt' dN ^\ M ) (35) 
NlK> J df S Nl (t',t) 1 ; 

ti 

which includes the relativistic correction factor for the decays of Ni as function of time. 

In addition to the baryon asymmetry we are interested in thermal production of grav- 
itinos which is governed by the Boltzmann equation 

d 

da G G 

with the dominant QCD collision term in the supersymmetric Standard Model [7|8] 



aH—N d = a A C d {T) , (36) 



mj \ 54C(3)^(T) 2 
L&{1 ] ~ 1 1 + 3m| J vr^Mp J 



2 



ln( ^^yr 1+0.8846 



(37) 



here ^(lO 10 GeV) = 0.85 is the QCD gauge coupling, and m g (T) = y^3/2g s (T)T is the 
plasma mass of the gluon. 

Finally, we need the temperature as function of the scale factor in order to compute 
the collision term for gravitino production and the equilibrium number density of the 
N\ neutrinos. The covariant energy conservation connects the thermal contributions to 
energy density and pressure with the nonthermal contributions, 

j t (pR + pl 1 +Ps + p S Nl )+3H(p R + p T Ni + Ps + p s Ni +p R + p T Ni +p s Ni )=0 . (38) 

The quantities p$, p Xl and p Xi are easily obtained from the corresponding distribution 
functions; energy and pressure of the thermally produced N\ neutrinos are given in terms 
of the number density, 

aY Nl =NlT, a *p T Ni =T^N T Ni . (39) 

Together with the Boltzmann equations (|26|) , (|34|) and (|36|) we now have a complete system 
of first-order differential equations which determine the time evolution of the initial state, 
in particular the generation of baryon asymmetry and gravitino abundance. 
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4 Results and discussion 



We have numerically solved the Boltzmann equations derived in the previous section. In 
the following, we describe the solution for a representative choice of neutrino masses, CP 
asymmetries^ and superparticle masses: ffi\ = 10~ 3 eV, Mi = 10 10 GeV, e\ = 10 -6 , 
£2,3 = —3 x 10~ 4 , and = 100 GeV, m g = 800 GeV. Technical details of the solution 
and a systematic study of the parameter space compatible with leptogenesis and gravitino 
dark matter will be presented in a forthcoming publication [H] . 

The various components of the energy density and the comoving number densities are 
compared as functions of the scale factor a in Fig. 1 and Fig. 2, respectively. Immediately 
after the end of tachyonic preheating the decay of N2 and N$ generates radiation, an initial 
B — L asymmetry and subsequently a thermal abundance of N\ neutrinos. At a ~ 2, the 
S bosons decay. The produced N\ neutrinos are relativistic, but nonthermal. At a ~ 150, 
the thermal Ni neutrinos reach thermal equilibrium, the initial B— L asymmetry is washed 
out, a new B — L asymmetry and the dominant part of radiation are produced in decays 
of the nonthermal Nf neutrinos. Furthermore, gravitinos are continuously produced from 
the thermal bath. Around a ~ 1000 the final B — L asymmetry and the gravitino number 
density are reached. 

The present value of baryon asymmetry is obtained from 



n° B g°N B - L 
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(40) 



where a/ ~ 10 4 is a scale factor after leptogenesis is completed; in the supersymmetric 
Standard Model the sphaleron conversion factor c sp h = 8/23, the effective numbers of 
degrees of freedom at high and low temperatures are g* = 915/4 and g® = 43/11, respec- 
tively, and the number density of photons is iV 7 = a 3 g 7 C(3)/7r 2 T 3 . The corresponding 
expression for the present gravitino abundance is given by 



n 



7 



(41) 



where p° c = 1.052 x 10 5 h? GeVcm 3 is the critical density. For our choice of parameters 
we obtain 

7? B = 1.6 x KT 7 , fi 6 /i 2 = 0.11. (42) 



6 We have chosen opposite signs for the CP asymmetries t\ and £2,3, so that one can distinguish their 
contribution to the final B — L asymmetry. 
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Figure 1: Comoving energy densities for S bosons, radiation (R), N% neutrinos produced in S decays 
(Nf), N\ neutrinos in equilibrium (N° q , for comparison), thermally produced Ni neutrinos (Nf) and 
gravitinos (G) as functions of the scale factor a. 




1 10 100 1000 io 4 



Scale factor a 



Figure 2: Comoving number densities for S bosons, B — L charge (B ~ L) N\ neutrinos produced in S 
decays (Nf), Ni neutrinos in equilibrium (N^ q , for comparison), thermally produced Ni neutrinos (Nj) 
and gravitinos (G) as functions of the scale factor a. 
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Figure 3: Radiation temperature as function of the scale factor a. 

The gravitino abundance is close to the observed value for DM . The calculated baryon 
asymmetry is about two orders of magnitude larger than the observed one, r]g MB ~ 
6.2 x 10~ 10 [10], which is consistent since we have used the maximal CP asymmetry. 

The temperature of radiation is shown in Fig. 3 as function of the scale factor. It 
decreases like T oc l/a except for an approximate plateau around Tl — 5 x 10 9 GeV 
at a ~ 100, when the N\ neutrinos decay. Note that the temperature Tl agrees with 
the reheating temperature Tr estimated in Eq. (j3J) up to a factor of two. Inserting 
the temperature Tl in the estimate Eq. (JTJ for the gravitino abundance, one obtains 
agreement with the result of the full calculation given in Eq. ( f42|) . For such large reheating 
temperatures nonthermal production of gravitinos is usually negligible |21j . 

It is instructive to compare the calculated baryon asymmetry with two extreme cases, 
namely thermal leptogenesis and the rapid conversion of a gas of nonrelativistic N\ neutri- 
nos dominating the energy density of the universe. For thermal leptogenesis one obtains 
for our choice of parameters (cf. [TT] ) 

^thermal = /( ^) ~ 5 x IO" 10 , (43) 

9* 

where we have used c sp h = 8/23 and kj(10~ 3 eV) ~ 0.1 for the efficiency factor. In the 
case of rapid conversion, energy conservation yields (cf. |22j) 

% Pid ^7^ sphei ^^9xlO- 7 , (44) 
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which is closer to the result of our calculation Eq. (142 p . Varying the parameters rh\ 
and Mi one can interpolate between thermal and nonthermal leptogenesis. The effective 
neutrino mass rh\ is closely related to the smallest neutrino mass mi and therefore to the 
absolute neutrino mass scale. Via its effect on the reheating temperature and gravitino 
dark matter, the absolute neutrino mass scale is thus related to the gravitino mass. 

The case of rapid conversion is also realized in leptogenesis from inflaton decays [23] . 
Here the reheating temperature is determined by the decay width of the inflaton. One 
can then have Tr <C Mi and in this way avoid overproduction of gravitinos. For larger 
inflaton decay widths reheating temperatures Tr ~ Mi can be reached, associated with 
larger baryon asymmetries [23], similar to the situation described in this paper. 

We have shown that reheating the universe by heavy Majorana neutrino decays can 
simultaneously explain the cosmological baryon asymmetry and the observed dark matter 
abundance in terms of thermally produced gravitinos. Starting from an initial state of 
unbroken B — L symmetry, tachyonic preheating leads to an interplay of nonthermal 
and thermal leptogenesis, which is controlled by the heavy neutrino decay widths. Open 
questions concern the connection with models of inflation and tachyonic preheating in 
supersymmetric theories. 
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